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Abstract. The measure and conquer approach has proven to be a powerful tool to analyse 
exact algorithms for combinatorial problems, like Dominating Set and Independent 
Set. In this paper, we propose to use measure and conquer also as a tool in the design of 
algorithms. In an iterative process, we can obtain a series of branch and reduce algorithms. 
A mathematical analysis of an algorithm in the series with measure and conquer results in 
a quasiconvex programming problem. The solution by computer to this problem not only 
gives a bound on the running time, but also can give a new reduction rule, thus giving 
a new, possibly faster algorithm. This makes design by measure and conquer a form of 
computer aided algorithm design. 

When we apply the methodology to a Set Cover modelling of the Dominating Set 
problem, we obtain the currently fastest known exact algorithms for Dominating Set: 
an algorithm that uses 0(1.5134") time and polynomial space, and an algorithm that uses 
0(1.5063") time. 



1. Introduction 

The design of fast exponential time algorithms for finding exact solutions to NP-hard 
problems such as Independent Set and Dominating Set has been a topic for research 
for over 30 years, see e.g., the results on Independent Set in the 1970s by Tarjan and 
Trojanowski [18, 19]. A number of different techniques have been used for these and other 
exponential time algorithms [5, 21, 22]. 

An important paradigm for the design of exact algorithms is branch and reduce, pi- 
oneered in 1960 by Davis and Putnam [1]. Typically, in a branch and reduce algorithm, 
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a collection of reduction rules and branching rules arc given. Each reduction rule simpli- 
fies the instance to an equivalent, simpler instance. If no rule applies, the branching rules 
generate a collection of two or more instances, on which the algorithm recurses. 

An important recent development in the analysis of branch and reduce algorithms is 
measure and conquer, which has been introduced by Fomin, Grandoni and Kratsch [4]. The 
measure and conquer approach helps to obtain good upper bounds on the running time of 
branch and reduce algorithms, often improving upon the currently best known bounds for 
exact algorithms. It has been used successfully on Dominating Set [4], Independent Set 
[6], Dominating Clique [15], the number of minimal dominating sets [9], Connected 
Dominating Set [7], Minimum Independent Dominating Set [12], and possibly others. 

In this paper, we show that the measure and conquer approach can not only be used for 
the analysis of exact algorithms, but also for the design of such algorithms. More specifically, 
measure and conquer uses a non-standard size measure for instances. This measure is based 
on weight vectors, which are computed by solving a quasiconvex programming problem. 
Analysis of the solution of this quasiconvex program yields not only an upper bound to the 
running time of the algorithm, but also shows where we should improve the algorithm. This 
can lead to a new rule, which we add to the branch and reduce algorithm. 

We apply this design by measure and conquer methodology to a Set Cover modelling 
of the Dominating Set problem, identical to the setting in which measure and conquer 
was first introduced. If we start with the trivial algorithm, then, we can obtain in a number 
of steps the original algorithm of Fomin et al. [4], but with additional steps, we obtain a 
faster algorithm, using 0(1.5134'*) time and polynomial space, with a variant that uses 
exponential memory and 0(1.5063"^) time. We also show that at this point we cannot 
improve this measure and conquer computed running time, unless we choose a different 
measure or change the branching rules. 

While for several classic combinatorial problems, the first non-trivial exact algorithms 
date many years ago, for the Dominating Set problem, the first algorithms with running 
time 0*(c") with c < 2 are from 2004, with three independent papers: by Fomin et al. [10], 
by Randerath and Schiermeyer [16], and by Grandoni [13]. The so far fastest algorithm for 
Dominating Set was found in 2005 by Fomin, Grandoni, and Kratsch [4]: this algorithm 
uses 0(1.5260") time and polynomial space, or 0(1.5137^*) time and exponential space. 

2. Preliminaries 

Given a collection of non-empty sets S, a set cover of (S is a subset C C. S such that 
every element in any of the sets in S occurs in some set in C. In the Set Cover problem 
we are given a collection S and are asked to compute a set cover of minimum cardinality. 

The universe Us of a Set Cover problem instance is the set of all elements in any set 
in <S; Us = UsgS ^- "^^^ frequency /(e) of an element e eUs is the number of sets in S in 
which this element occurs. Let >S(e) = {S G >S | e G S"} be the set of sets in S in which the 
element e occurs. Wc define a connected component C in a Set Cover problem instance 
<S in a natural way: a minimal non-empty subset C C. S for which all elements in the sets in 
C occur only in sets contained in C. The dimension ds of a Set Cover problem instance 
is the sum of the number of sets in S and the number of elements in Us', ds = \S\ + \Us\- 

Let G = (y, E) be an undirected graph. A subset D CV of nodes is called a dominating 
set if every node v E V is either in D or adjacent to some node in D. The Dominating 
Set problem is to compute for a given graph G a dominating set of minimum cardinality. 
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We can reduce the minimum dominating set problem to the Set Cover problem by 
introducing a set for each node of G which contains the node itself and its neighbours; 
S := {-/V[i'] I V G y}. Thus we can solve a Dominating Set problem on a graph of n nodes 
by a minimum set cover algorithm running on an instance of dimension d = 2n. 

Both problems are long known to be NP-complete [11, 14], which motivates the search 
for fast exponential time algorithms. 

3. A Faster Exact Algorithm for Dominating Set 

In this section, we give our new algorithm for Dominating Set. The algorithm is 
an improvement to the algorithm by Fomin et al. [4]; it is obtained from this algorithm 
by adding some additional reduction rules. These rules were derived using the design by 
measure and conquer approach, see Section 4. After introducing our algorithm, we recall 
the necessary background of the measure and conquer method [4] . 

3.1. The Algorithm 

Our algorithm for the DOMINATING Set problem uses the Set Cover modelling of 
Dominating Set shown in Section 2. It is a branch and reduce algorithm on this modelling 
consisting of four simple reduction rules, one base case for the recursion and a branching 
rule. See Algorithm 1. 

For a given problem instance we first apply the following reduction rules: 

(1) Base case. If all sets in the instance arc of size at most two then finding a minimum 
set cover is equivalent to finding a maximum matching in a graph. Introduce a node 
for each element and an edge for each set of size two. Now the maximum matching 
joined with some sets containing the unmatched nodes form a minimum set cover. 
This matching can be computed in polynomial time [2]. 

(2) Splitting connected components. If the instance contains multiple connected compo- 
nents, compute the minimum set cover in each connected component separately. 

(3) Subset rule. If the instance contains sets Si, S2 with Si C 5*2, then we remove 5*1 
from the instance. Namely, in each minimum set cover that contains ^i, we can 
replace by S2 and obtain a minimum set cover without S2. 

(4) Subsumption rule. If the set of sets ^(e') in which an element e' occurs is a subset 
of the set of sets >S'(e) in which another element e occurs, we remove the element e. 
For any set cover, covering e' also covers e. 

(5) Unique element or singleton rule. If any set of size one remains in the instance after 
application of the previous rules, we add this set to the set cover. For the element 
in this set must occur uniquely in this set, otherwise it would have been a subset of 
another set and have been removed by rule 3. 

(6) Avoiding unnecessary branchings based on frequency two elements. For any set S in 
the problem instance let r2 be the number of frequency two elements it contains. Let 
m be the number of elements in the union of sets containing the other occurrences 
of these frequency two elements, excluding any element already in S. If for any set 
S: m < r2 then include S in the set cover. 

This rule is correct since if we would branch on S and include it in the set cover we 
would cover \S\ elements with one set. If the set cover does not use S, it must use all 
sets containing the other occurrence of the frequency two elements, since they have 
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become unique elements now. Notice that by Rule 4 all other occurrences of the 
frequency two elements must be in different sets and thus we would cover \S\ + m 
elements with r2 sets. So if m < r2 the first case can be preferred over the second: 
we can just add S to the cover and have r2 — 1 < m sets left to cover at least this 
much elements. 

For the branching rule, we select a set of maximum cardinality and create two subproblems 
by either including it in the minimum set cover and removing all newly covered elements 
from the problem instance or removing it. 

Algorithm 1 Algorithm Designed by Measure and Conquer 
MSC(5) = { 

if max{\S\ \S eS}<2 then 

return minimum set cover of S by computing a matching 
else if 3C C 5 : C ^ {0, 5}, {S{e)\e eS,S eC} = C then 

return MSC(C) + MSC(5\C) 
else if 3S,S' eS : S ^ S', S C S' then 

return MSC(5\{5}) 
else if 3e, e' eUs-e^ e', <S(e) D S{e') then 

S' = MSC({5\{e}|5 G S}) 

return {S U {e}\S G 5', 5 U {e} G 5} U {S\S G 5', 5 U {e} ^ <S} 
else if 3{e} G S then 

return {{e}} U MSC(5\{{e}}) 
else if 3S : |U{5'\'5|e G S,f{e) = 2,S' E S{e)}\ < \{e G S\f{e) = 2}| then 

return S U MSC({5'\5|5' G S\{S}}) 
else 

Let S := argmax^/g^dS"!) 

r = {MSC(5\{5}), S U MSC({5'\5|5' G 

return argminpgp(|P|) 
end if 
} 



3.2. Running time analysis by Measure and Conquer 

The basic idea of measure and conquer is the usage of a non-standard measure for the 
complexity of a problem instance in combination with an extensive subcase analysis. In the 
case of Set Cover, we give weights to set sizes and element frequencies, and sum these 
weights over all items and sets. We enumerate many subcases in which the algorithm can 
branch and derive recurrence relations for each of these cases in terms of these weights. 
Finally we obtain a large numerical optimisation problem which computes the weights 
corresponding to the smallest solution to all recurrence relations, giving an upper bound on 
the running time of our algorithm. This analysis is similar to [4] . 

We let Vi,Wi G [0,1] be the weight of an element of frequency i and a set of size i 
respectively, and set our variable measure of complexity ks to: 

ks = ^ w\s\ + ^ w/(e) notice : ks < ds 
ses eeUs 
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Sets of different sizes and elements of different frequencies contribute equally to the dimen- 
sion of the instance, but now larger sets and higher frequency elements can contribute more 
to the measured complexity of the instance. Furthermore we set Vi = Wi = 0,i G {0,1} 
since all frequency one elements and size one sets are removed by the reduction rules. For 
later use we introduce quantities representing the reduction in problem complexity when 
the size of a set or the frequency of an element is reduced by one. For technical reasons, 
these quantities must be non-negative. 

Awi = Wi — Wi-i Avi = Vi — Vi-i Vz > 1 : Avi, Awi > 

The next step will be to derive recurrence relations representing problem instances the 
algorithm branches on. Let N{k) be the number of subproblems generated in order to solve 
a problem of measured complexity k. And let Akin (include S) and Akout (discard S) 
be the difference in measured complexity of both subproblems compared to the problem 
instance we branch on. Finally let l^l = Yli^2 where the number of elements in S of 
frequency i. 

If we add S to the set cover, S is removed together with all its elements. This re- 
sults in a reduction in size of w\s\ + J2i^2''^i'^i- Because of the removal of these ele- 
ments, other sets are reduced in size; this leads to an additional complexity reduction 
of at least mm j^^g^{ Aw j}'^^2i'^ ~ To keep the formula (and the next) linear, we 

set minj<|5|{Ait;j} = Aw^s] ^^'^ keep the formula correctly modelling the algorithm by 
introducing the following constraints on the weights: 

Vz > 2 : Awi < Awi-i 

One can show that including these in the numerical optimisation problem does not change 
the solution, as it gives the same weights. The constraints help to considerably speed up 
this optimisation process. 

If we discard S, we also remove it from the problem instance, and hence all its elements 
are reduced in frequency by one. So we have a complexity reduction of w\s\ + X]£2 ^i^'^^i- 
Besides this reduction, the sets which contain the second occurrences of any frequency 
two element are included in the set cover. Notice that these must be different sets due to 
reduction Rule 4. Because of Rule 6 we know that at least r2 other elements must be in these 
sets as well, and these also must occur somewhere else in the instance, hence even more sets 
are reduced in size. Summation leads to an additional size reduction of r2{v2 + W2 + Aw^s\)- 
Here we also use Rule 2 to make sure that not all these frequency two elements occur in the 
same set, because in that case all considered sets form a connected component of at most 
five sets which thus can be solved in 0(1) time. 

This leads to the following set of recurrence relations: V 3 < l^*! = X^^2 

N{k) < N{k - Akout) + N{k - Akin) 

where 

CO 

Afcowt = w^si + ^ nAvi + r2{v2 + W2 + Aw^s\) 

i=2 

oo oo 

Akin = W^|5| + XI + ^'^\S\ ^i'i' - ^)^i 
1=2 i=2 

We make the problem finite by setting for some large enough p all vi = Wi = 1 for 
i > p, and only consider the subcases l^l = Yli=2''^i + ''>P' where r>p = Yli^p+i''^i- ^ow 
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we have a finite set of recurrences which model our algorithm since the recurrences for the 
cases where \S\ > p + 1 arc dominated by those where \S\ = p + I. The best value for p 
follows from the optimisation, for if chosen too small the now constant recurrences (weights 
equal 1) will dominate all others in the optimum, and if chosen too large the extra Vi and Wi 
are optimised to 1 (and the optimisation problem was unnecessarily hard). Here p equals 7. 

A solution to this set of recurrence relations will be of the form N{k) = , where a is 
the smallest solution of the set of inequalities: 

Since k < d where d the dimension of the problem, we know that the algorithm will have a 
running time of 0{{a + e^), for any e > 0: 

0{poly{d)N{k)) = 0{poly{d)Q^) < 0{poly{d)a'^) < 0((a + e)'^) 

Prom here on we let e be the error in the upward decimal rounding of a. 

So for any given vector v = (0, V2. v^.^ V/i, . . .) and w = (0, W2,w^, W/i, . . .) we can now 
compute the running time measured with these weights. As a result we have obtained a 
numerical optimisation problem: choose the best weights so that the upper bound on the 
running time is minimal. 

The numerical solution to this problem can be found in the last cell of Table 1, resulting 
in an upper bound on the running time of the algorithm of 0(1.2302''): 

3.3. Quasiconvex programming 

The sort of numerical optimisation problems arising from measure and conquer analyses 
are quasiconvex programs, named after the kind of function we are optimising: a quasiconvex 
function, which is a function with convex level sets {x \ q{x) < A}. 

To our knowledge there are currently two different techniques in use to solve these 

quasiconvex programs: randomised search, and Eppstcin's smooth quasiconvex programming 
algorithm [3]. We have implemented a variant of the second technique; for details see [20]. 

3.4. Results 

As discussed, we have now obtained the following result. 

Theorem 3.1. Algorithm 1 solves a Set Cover problem instance of dimension d in 
0(1.2302*^) time and polynomial space. m 

Using the minimum set cover modelling of Dominating Set this results in: 

Corollary 3.2. There exists an algorithm that solves the Dominating Set problem in 
0(1.5134") time and polynomial space. m 

We can further reduce the time complexity of the algorithm at the cost of exponential 
space. This can be done by dynamic programming; the algorithm keeps track of all solutions 
to subproblems solved and if the same subproblem turns up more than ones it is looked up. 
Notice that querying and storing the subproblems can be implemented in polynomial time. 

We compute the new time complexity based on [8, 17] and obtain: 



Theorem 3.3. Algorithm 1, modified as above, solves a Set Cover problem of dimension 
d in 0(1.2273"^) time and space. ■ 
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Corollary 3.4. There exists an algorithm that solves the DOMINATING Set problem in 



4. Design by Measure and Conquer 

The beauty of our algorithm hes in the fact that it has been designed using a form 
of computer aided algorithm design which we call design by measure and conquer. Given 
a variable measure of complexity as in the analysis in Section 3.2 and a set of branching 
rules, all polynomial time computable reduction rules relative to this measure and branching 
rules follow by the method. We start with a trivial branch and reduce algorithm, i.e. one 
without any reduction rules and only consisting of the branching rule and a trivial base case 
(if the problem is empty, return 0). Next we exhaustively apply an improvement step, which 
comes up with a new reduction rule and hence a possibly faster algorithm. This changes the 
algorithm analysis technique measure and conquer into a technique for algorithm design. 

Thus, this gives a very nice process, where a human invents additional reduction rules, 
and the computational power of our computer does the extensive measure and conquer 
analysis and points to all possible points of direct improvement. This combination has 
proven to be successful as we see from the results of Section 3. While constructing our 
algorithm, the previously fastest algorithm for DOMINATING Set by Fomin et al. [4] has 
been obtained as an intermediate step. It has now been improved up to a point where 
we need to either change the branching rule (or add new branching rules) or modify the 
measure and conquer analysis, i.e. use a different variable measure or perform a more 
elaborate subcase analysis. See Table 1 for information on the analysis and added rules for 
all algorithms, from the starting trivial algorithm without any reduction rules till we obtain 
Algorithm 1. 

4.1. A Single Iteration: improving the previously fastest algorithm 

We now demonstrate how the improvement step works, by giving one such improve- 
ment as elaborate example, namely an improvement we can make when we start with the 
algorithm by Fomin ct al. [4]. This step is marked with a star in Table 1. First wc perform a 
measure and conquer analysis on the current algorithm giving us the optimal instantiation 
of the variable measure, and an upper bound on the running time of the algorithm. Next 
we examine the quasiconvex function we have just optimised. 

The quasiconvex function has the following form: 



where C is the set of all possible instances the algorithm can branch on and AA;^„, A/c„,^j^ arc 
the differences in measured complexity between the generated subproblems and branching 
subcase c. 

Each one of the functions Qc is quasiconvex (see [3]), i.e. it has convex level sets. The 
situation is very similar to finding the point x of minimum maximum distance to a set of 
points P in N dimensional space: only a few points in P have distance to x tight to this 
maximum, and moving away from x always results in at least one of these distances to 
increase. If one such tight point is moved or removed, this directly influences the optimum 
X and the minimum maximum distance. 



0(1.5063") time and space. 
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We now consider the eight cases that are tight to the value of the quasiconvex function 
q in the optimum. These are: 

|5|=r-2 = 3, |5| = r3 = 3, 1^1 = r4 = 3, \S\ = = A 
1^1= re = 4, |5|=r6 = 5, \S\ = rj = f,, \S\ = tt = Q 

Now, if we can formulate a reduction rule that either further reduces the size of any sub- 
problem generated in these cases, or removes any of these cases completely, then we lower 
the value of the corresponding Qc, or remove this Qc respectively, resulting in a new optimum 
corresponding to a faster running time. 

We take the simplest case for improvement; IS"! as small as possible, and with as low 
frequency elements as possible. This corresponds to an instance with: 

5 = {1,2, 3} existing next to: { 1, 2, 4 }, { 3, 4 } 

We emphasize that this is not an entire instance, but just a fragment of an instance con- 
taining the set S used for branching and the collection of sets in which the elements from S 
also occur. In an instance corresponding to this subcase the element 4 can be of frequency 
two or higher, but all sets are of size three or smaller. 

We note that we do not need to branch on this particular subcase: elements 1 and 2 
occur in exactly the same sets, and thus if a set cover covers one of these, the other is covered 
as well. We generalise this and formulate the subsumption rule (Rule 4 of Algorithm 1). 
Now we have a new algorithm, for which we can adjust the measure and conquer analysis, 
and repeat this process. 

4.2. The Process Halts 

Above, we discussed how to perform one step of the design by measure and conquer 
process. For a complete overview of the construction of Algorithm 1 see Table 1, with the 
relevant data for each improvement step. Note from Table 1 that after each new step, the 
example worst case instance part is no longer a valid worst case for the next step. As a 
result, at each step either some subcases are removed by using a larger smallest set S or 
by removing small sets or elements (setting = or = 0), or the size reduction in 
the formula for Akout is increased. After each step we refactored the reduction rules and 
removed possible superfluous ones. We have not included the formula for Akin in this table, 
since it does not change except that ri 7^ in early stages. 

It appears that we must use a different approach to obtain a faster algorithm. Consid- 
ered the following problem: 

Problem 4.1. Given a Set Cover instance S and a set 5 € <S with the properties: 

(1) Non of the reduction rules of Algorithm 1 apply to >S. 

(2) All sets in S have cardinality at most three; |5| = 3. 

(3) Every element e e S has frequency two. 

Question: Does there exist a minimum set cover of S containing S? 

Proposition 4.2. Problem 4-1 is NP-complete. 

Proposition 4.2 implies that we cannot formulate a polynomial time reduction rule that 
removes the current simplest worst case of our algorithm by deciding on whether S" is in a 
minimum set cover or not, unless P = NP. 
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Latest new reduction rule 

current formula for Akout 
subcases considered 
weights vectors v and w 



Running times for Set Cover and 

Dominating Set 
instance part of the simplest worst case; 

S — other occurrences of elements of S 



Trivial algorithm 

l<\S\=T,i=iri + r>p<p+l = 3 
v = (0.8808,0.9901, . . .) 



0(1.4519'*) 0(2.1080") 

{1} - 

w = (0.9782,...) 



Stop when all sets of size one 

w\s\ + E£i n^vi 

2<\S\ = Ef=i + r>p < p + 1 = 4 
v = (0.7289,0.9638,0.9964,...) 



0(1.3380'*) 0(1.7902" 

{1,2} - 

w = (0.4615,0.9229,...; 



Include all frequency one elements 0(1.2978'*) 0(1.6842") 

{1,2} -{1,2} 

w = (0.4240, 0.8480, 0.9676, . . .) 



2<\S\ = Zt^ri + ryp<p+l = 5 
V = (0, 0.4818, 0.8357, 0.9636, . . .) 



Subset rule 



0(1.2665'*) 0(1.6038") 



+ + '^r2>o(w2 + ^2) + 5\s\=r2=2^W2 

, . ^..i-'z {1,2} - {1,3},{1,4},{2,3},{2,4} 

w = (0, 0.6973, 0.9093, 0.9800, . . .) 



2< |5| = Er=2^« + ^>P<P + l = 7 
V = (0, 0.3900, 0.7992, 0.9318, 0.9808, . . . ) 



Compute matching for size two sets* 0(1.2352^ 0(1.5258") 

+ Ei^^i^^i + '^r2>o('i^2 + '"2) + 5^s\=i,r2>2{^Wz + ^^2=3^2) + 5^s\=r2=i'^i 

{1,2,3} - {1,2,4}, {3,4} 
w = (0, 0.3787, 0.7575, 0.9103, 0.9763, 



3< |5| = E"=2^i + ^>P<P + l = 7 
V = (0, 0.3978, 0.7650, 0.9263, 0.9842, . . .) 



0(1.2339**) 0(1.5223") 



Subsumption rule 



W\S\ + ES2 ^i^Vi + 5r2>Q{r2W2 + V2) + (5|5|=^2=3 A?;3 

3<\S\ = j:t2n + r>p<p+l = 7 {1,2,3} - {1,4},{2,4},{3,4} 

V = (0, 0.3545, 0.7455, 0.9203, 0.9818, . . .) w = {0, 0.3755, 0.7510, 0.9061, 0.9745, 
Avoid unnecessary branchings 0(1.2313**) 0(1.5160") 

+ Ei^2''i^^i + ^^2(^2 + V2) + Sr2>l{r2 - l)Au'|s| 

3<\S\ = Et2ri + r>p<p+l = 8 {1,2,3} - {1,4},{2,5},{3,6} 

V = (0, 0.269912, 0.689810, 0.892666, 0.965849, 0.992140, . . .) 
w = (0, 0.376088, 0.752176, 0.907558, 0.974394, 0.999212, . . .) 



Connected components (final) O(1.2302<*) 0(1.5134") 

W\S\ + Y^'^2'^i^^i + ^2(^2 +V2 + Awis\) 

3<\S\ = TH=2ri + r>p<p+l = 8 {1,2,3} - { 1, 4 }, { 2, 5 }, { 3, 6 } 

V = (0, 0.219478, 0.671386, 0.876555, 0.956850, 0.988195, . . .) 

w = (0, 0.375418, 0.750835, 0.905768, 0.971965, 0.998158, . . .) 

* Algorithm by Fomin, Grandoni and Kratsch [4]. 

Table 1: The iterations of the design by measure and conquer process. 
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We can construct similar NP-complete problems for all other worst cases of Algorithm 1. 
Therefore Algorithm 1 is optimal in some sense: wc cannot straightforwardly improve it by 
performing another iteration. In order to obtain a faster branch and reduce algorithm using 
polynomial time reduction rules with a smaller measure and conquer proved time bound, 
it is necessary to either change the variable measure, the branching rule(s), or perform a 
more extensive subcase analysis. 

Very recently, we pursued the last option with little result. We tried to further subdivide 
the frequency two elements in the branch set depending on the size of the set containing 
their second occurence (two or larger) and if this is a set of size two, on the frequency of 
the other element in this set. This resulted in a set of very technical reduction rules and a 
small speedup for the case where we use only polynomial space. This speedup, however, was 
almost completely lost when using exponential space because some of the weights involved 
became almost zero. 

5. Conclusion and Further Research 

In this paper, we have given the currently fastest exact algorithm for the Dominating 
Set problem. Besides setting the current record for this central graph theoretic problem, 
we also have shown that measure and conquer can be used as a tool for the design of 
algorithms. 

We have shown that there exists a strong relation between the chosen variable measure, 
the branching rulc(s) and the reduction rules of a measure and conquer based algorithm. 
We intend to further investigate this relation and examine to what point we can deduce not 
only reduction rules, but also branching rules from the given measure. 

We plan to apply the design by measure and conquer method to a number of other 
combinatorial problems, and hope and expect that in a number of cases, such a computer 
aided algorithm design will give further improvements to the best known exact algorithms 
for these problems. 

In this paper, we observe that measure and conquer can be used as a form of computer 
aided algorithm design. Another intriguing question is whether we can automate some 
additional steps in the design process, e.g., can we automatically obtain reduction rules 
from the solution of the quasiconvex program? 
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